Single-shot measurement of quantum optical phase 
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Although the Canonical phase of light, which is defined as the complement of photon number, 
has been described theoretically by a variety of distinct approaches, there have been no methods 
proposed for its measurement. Indeed doubts have been expressed about whether or not it is mea- 
surable. Here we show how it is possible, at least in principle, to perform a single-shot measurement 
of Canonical phase using beam splitters, mirrors, phase shifters and photodetectors. 

PACS numbers: 42.50.Dv, 42.50.-p 
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Quantum-limited phase measurements of the optical 
field have important applications in precision measure- 
ments of small distances in interferometry and in the 
emerging field of quantum communication, where there 
is the possibility of encoding information in the phase of 
light pulses. Much work has been done in attempting 
to understand the quantum nature of phase. Some ap- 
proaches have been motivated by the aim of expressing 
phase as the complement of photon number []J . Exam- 
ples of these approaches include the probability operator 
measure approach]^, a formalism in which the Hilbert 
space is doubled [4, a limiting approach based on a fi- 
nite Hilbert space ja, EJ and a more general axiomatic ap- 
proach 0. Although these approaches are quite distinct, 
they all lead to the same phase probability distribution 
for a field in state \ip) as a function of the phase angle 9 
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where \n) is a photon number state. Leonhardt et al. 
have called this common distribution the "canoni- 
cal" phase distribution to indicate a quantity that is 
the canonical conjugate, or complement, of photon num- 
ber. This distribution is shifted uniformly when a phase 
shifter is applied to the field and is not changed by a 
photon number shift lj. We adopt this definition here 
and use the term Canonical phase to denote the quantity 
whose distribution is given by Q). 

Much less progress has been made on ways to measure 
Canonical phase. Homodyne techniques can be used to 
measure phase-like properties of light but are not mea- 
surements of Canonical phase. It is possible in principle 
to measure the Canonical phase distribution by a series of 
experiments on a reproducible state of light but there 
has been no known way of performing a single-shot mea- 
surement. Indeed it is thought that this might be impos- 
sible [U- Even leaving aside the practical issues, the con- 
cept that a particular fundamental quantum observable 
may not be measurable, even in principle, has interesting 
general conceptual ramifications for quantum mechanics. 
A different approach to the phase problem, which avoids 



difficulties in finding a way to measure Canonical phase, 
is to define phase operationally in terms of observables 
that can be measured . The best known of these op- 
erational phase approaches is that of Noh et al. [3, . 
Although the experiments to measure this operational 
phase produce excellent results, they were not designed 
to measure Canonical phase as defined here and, as shown 
by the the measured phase distribution 0, they do not 
measure Canonical phase. In this paper we show how, 
despite these past difficulties, it is indeed possible, at 
least in principle, to perform a single-shot measurement 
of Canonical phase in the same sense that the experi- 
ments of Noh et al. are single-shot measurements of their 
operational phase. 

A single-shot measurement of a quantum observable 
must not only yield one of the eigenvalues of the ob- 
servable, but repeating the measurement many times on 
systems in identical states should result in a probability 
distribution appropriate to that state. If the spectrum 
of eigenvalues is discrete, the probabilities of the results 
can be easily obtained from the experimental statistics. 
Where the spectrum is continuous, the probability den- 
sity is obtainable by dividing the eigenvalue range into a 
number of small bins and finding the number of results 
in each bin. As the number of experiments needed to ob- 
tain measurable probabilities increases as the reciprocal 
of the bin size, a practical experiment will require a non- 
zero bin size and will produce a histogram rather than a 
smooth curve. 

Although the experiments of Noh et al. 0, E| were 
not designed to measure Canonical phase, it is helpful 
to be guided by their approach. In addition to their re- 
sults being measured and plotted as a histogram, some 
of the experimental data are discarded, specifically pho- 
ton count outcomes that lead to an indeterminacy of the 
type zero divided by zero in their definitions of the cosine 
and sine of the phase Q . The particular experiment that 
yields such an outcome is ignored and its results are not 
included in the statistics. 

Concerning the discarding of some results, we note in 
general that the well-known expression for the probabil- 
ity that a von Neumann measurement on a pure state \tp) 
yields a result q is (tp\q) (q\ip) , where \q) is the eigenstate 
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of the measurement operator corresponding to eigenvalue 
q. If the state to be measured is a mixed state with 
a density operator p the expression becomes Tr(/5II g ) 
where U q = \q)(q\. The operator II 9 is a particular case 
of an element of a probability operator measure (POM) 
0. The sum of all the elements of a POM is the unit 
operator and the expression Tr(/5II g ) for the probabil- 
ity is based on the premise that all possible outcomes 
of the measurement are retained for the statistics. If 
some of the possible outcomes of an experiment are dis- 
carded, the probability of a particular result calculated 
from the final statistics is given by the normalized ex- 
pression Tr(/5II 9 )/ 5Z p Tr(pIIp), where the sum is over all 
the elements of the POM corresponding to outcomes of 
the measurement that are retained. 

We seek now to approximate the continuous distribu- 
tion Q by a histogram representing the probability dis- 
tribution for a discrete observable 9 m such that when the 
separation 59 of consecutive values of 8 m tends to zero 
the continuous distribution is regained. A way to do this 
is first to define a state 
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There are N + 1 orthogonal states \6 m ) corresponding to 
N+l values m = mS6 with 59 = 2ir/(N + 1) and m = 
0, 1, ... N. This range for m ensures that 9 m takes values 
between and 27r. Then, if we can find a measurement 
technique that yields the result 9 m with a probability 
of KV'l^m)! 2 : the resulting histogram will approximate 
a continuous distribution with a probability density of 
I (V'l^m) \ 2 /S9. It follows that, as we let N tend to infinity, 
there will exist a value of 9 m as close as we like to any 
given value of 9 with a probability density approaching 
P{9) given by QJ. If we keep N finite so that we can 
perform an experiment with a finite number of outcomes, 
then the value of N must to be sufficiently large to give 
the resolution 59 of phase angle required and also for 
\ip) to be well approximated by X)n( n IV')l n ) where the 
sum is from n = to AT. The latter condition ensures 
that the terms with coefficients (n\ijj) for n > N have 
little effect on the probability |(^|# m }| 2 . As we shall be 
interested mainly in weak optical fields in the quantum 
regime with mean photon numbers of the order of unity, 
the maximum phase resolution 59 desired will usually be 
the determining factor in the choice of N. 

When N is finite, the states \9 m ) do not span the whole 
Hilbert space, so the projectors \9 m )(9 m \ will not sum to 
the unit operator 1. Thus these projectors by themselves 
do not form the elements of a POM. To complete the 
POM we need to include an element 1 — \8m)(9 m \. 
If we discard the outcome associated with this element, 
that is, treat an experiment with this outcome as an un- 
successful attempt at a measurement in a similar way 
that Noh et al. |9j treated experiments with indetermi- 
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FIG. 1: Multi-port device for measuring phase. The input 
and output modes are labelled 0,1, ... N from the top. In 
input mode is the field in state \tp)o to be measured and 
in input mode 1 is the reference field in state |-B)i. Vacuum 
states form the other inputs. There is a photodetector in each 
output mode. If all the photodetectors register one count 
except the detector D m in output mode m, which registers 
no counts, then the detector array acts as a digital pointer 
mechanism indicating a measured phase angle of 8 m . 



nate outcomes, then the probability that the outcome of 
a measurement is the phase angle 9 m is given by 
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where p = 0, 1 . . . N. We now require a single-shot mea- 
suring device that will reproduce this probability in re- 
peated experiments. 

As measurements will be performed ultimately by pho- 
todetectors, we seek an optical device that transforms 
input fields in such a way that photon number measure- 
ments at the output ports can be converted to phase 
measurements. We examine the case of a multi-port de- 
vice with N + l input modes and N + l output modes 
as depicted schematically in Fig. ^ As phase is not an 
absolute quantity, that is it can only be measured in re- 
lation to some reference state, we shall need the field in 
state \ip)a that is to be measured to be in one input and 
a reference field in state \B)i to be in another input. We 
let there be vacuum states |0)j with i = 2, 3 . . .N in the 
other input modes. We let the device be such that the 
input states are transformed by a unitary transforma- 
tion R. We let N + l photodetectors that can distinguish 
between zero photons, one photon and more than one 
photon be in the output modes. 

Consider the case where one photon is detected in each 
output mode except for mode m, in which zero photons 
are detected. The amplitude for this detection event is 
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where the first product is over values of j from to N 
excluding the value m. We require the transformation 
and the reference state to be such that this amplitude is 
proportional to o^mlV^o • We write the photon creation 
operators acting on the mode i as d\. Writing as 
o.^ 1 0) j- we see that we require \9 m )o to be proportional to 
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We rewrite the unitary transformation in the form 
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where Uja are the elements of a unitary matrix. Reck 
et al. 11] have shown how it is possible to construct a 
multi-port device from mirrors, beam splitters and phase 
shifters that will transform the input modes into the out- 
put modes in accord with any (TV + 1) x (N + 1) unitary 
matrix. We choose such a device for which the associated 
unitary matrix is 
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where u> = exp[— i2Tr/(N+ 1)] that is, a (iV + l)th root of 
unity. Such a device will transform the combined vacuum 
state to the combined vacuum state so we can delete the 
W on the right of expression JSJ. Substituting into 
© gives eventually 
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where k x = (N + l)- N / 2 . 

To evaluate (jSJ) we divide both sides of the identity 
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by X + co m to give, after some rearrangement and appli- 
cation of the relation oj m ( JV + 1 ) = 1, 



JJ(X + ^') = (-l) 



do) 



1 - (-Xcj- m ) 



The last factor is the sum of a geometric progression. Ex- 
panding this and substituting X = x/y gives eventually 
the identity 



N 



N- 



(11) 



71 = 



We now expand \B)± in terms of photon number states 



as 



\B) 1 = Y J K\n) 1 
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and put x = a J and y = a\ in (|ll|l . Then from JSJl we 
find that (JSJ becomes 

.<Bi(n«<°i) ( n^^i^ fni°)' 
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where k,2 = Kicj '"(iV!) 1 / 2 . We see then that, if we let 
\B)i be the binomial state 



iB) 1 = 2- w / 2 ^(-ir 



1/2 



(14) 



then expression (|13|l is proportional to X) n w ~™ m l n )oj 
that is, to |0 m )o- Thus the amplitude for the event that 
zero photons are detected in output mode m and one 
photon is detected in all the other output modes will 
be proportional to o(6 m \ip)o- The probability that the 
outcome of a measurement is this event, given that only 
outcomes associated with the (N +1) events of this type 
are recorded in the statistics, will be given by ©, where 
we note that the proportionality constant K2 will cancel 
from this expression. Thus the measurement event that 
zero photons are detected in output mode m and one 
photon is detected in all the other output modes can be 
taken as the event that the result of the measurement of 
the phase angle is 9 m . Thus the photodetector with zero 
photocounts, when all other photodetectors have regis- 
tered one photocount, can be regarded as a digital pointer 
to the value of the measured phase angle. 

We have shown, therefore, that it is indeed possible in 
principle to conduct a single-shot measurement of Canon- 
ical phase to within any given non-zero error, however 
small. This error is of the order 2tt/(N + 1) and will 
determine the value of N chosen. 

While the aim of this paper is to establish how Canoni- 
cal phase can be measured in principle, it is worth briefly 
considering some practical issues. Although we have 
specified that the photodetectors need only be capable 
of distinguishing among zero, one and more than one 
photons, reflecting the realistic case, there are other im- 
perfections such as inefficiency. These will give rise to 
errors in the phase measurement, just as they will cause 
errors in a single-shot photon number measurement. In 
practice, there is no point in choosing the phase reso- 
lution 66 much smaller than the expected error due to 
photodetector inefficiencies, thus there is nothing lost in 
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FIG. 2: Triangular array for N + 1 = 4. The outside beam 
splitters in the top row are 50:50, the middle one is fully 
reflecting. The second row beam splitters are |-transmitting 
and the bottom one is |-transmitting. The phase-shifter i 
produces a n/2 phase shift. 

practice in keeping N finite. A requirement for the mea- 
suring procedure is the availability of a binomial state. 
Such states have been studied for some time but their 
generation has not yet been achieved. In practice, how- 
ever, we are usually interested in measuring weak fields in 
the quantum regime with mean photon numbers around 
unity 

0,113 

and even substantially less 0|. Only the 
first few coefficients of \n)o in l|13[> will be important for 
such weak fields. Also, it is not difficult to show that 
the reference state need not be truncated at n = N, as 
indicated in 1)12(1 . as coefficients b n with n > N will not 
appear in {Tffi . Thus we need only prepare a reference 
state with a small number of its photon number state 
coefficients proportional to the appropriate binomial co- 
efficients. Additionally, of course, in a practical exper- 
iment we are forced to tolerate some inaccuracy due to 
photodetector errors, so it will not be necessary for the 
reference state coefficients to be exactly proportional to 
the corresponding binomial state coefficients. These fac- 
tors give some latitude in the preparation of the refer- 
ence state. The muti-port device depicted in Fig. 1 can 
be constructed in a variety of ways. Reck et al. 
provide an algorithm for constructing a triangular array 
of beam splitters to realize any unitary transformation 
matrix. Fig. 2 shows, for example, such a device that, 
with suitable phase shifters in the output modes, will 
realize the transformation Uij in JJJ) with N + 1 = 4. Be- 
cause we are detecting photons, however, these output 
phase shifters are not actually necessary and are thus 
not shown. The number of beam splitters needed for a 
general triangular array increases quadratically with N. 
Fortunately, however, our required matrix Q represents 
a discrete Fourier transformation and we only require two 
of the input ports to have input fields that are not in the 
vacuum state. The device of Torma and Jex is ide- 
ally suited for these specific requirements. This device, 
which has an even number of input and output ports, 



is pictured in Ref. [T^. It consists of just (N + l)/2 
ordinary 50:50 beam splitters and two plate beam split- 
ters. The latter are available in the form of glass plates 
with modulated transmittivity along the direction of the 
incoming beam propagation 14]. 

As a large fraction of the raw data in this procedure is 
discarded, an interesting question arises as to whether or 
not there is a relationship between the method of this pa- 
per and a limiting case of the operational phase measure- 
ments of Torgerson and Mandel where it is found that 
the distribution becomes sharper as more data are dis- 
carded. While preliminary analysis indicates that there 
is not, this will be discussed in more detail elsewhere. 

In conclusion we have shown that it is possible in prin- 
ciple to perform a single-shot measurement of the Canon- 
ical phase in the same sense that the experiments of Noh 
et al. are single-shot measurements of operational phase. 
The technique relies on generating a reference state with 
some number state coefficients proportional to those of a 
binomial state. 

D. T. P. thanks the Australian Research Council for 
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